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摘 要 :为 探讨 鱼 类 养殖 网 箱 对 周围 海域 流 场 的 影响 ,基于 雷诺 时 均 方 程 和 Realizable k-e 湛 流 模 
型 ,看 合 多 孔 介质 模型 ,建立 了 研究 鱼 类 养殖 网 和 和 周转 流 场 的 三 维 水 动力 模型。 根据 典型 养殖 场 的 
实际 数据 ,分 别 对 仅 网 衣 作 用 、 仅 和 鱼 类 作用 及 网 衣 和 和 鱼 类 共同 作用 下 网 箱 周 围 流 场 进行 数值 模拟 。 
结果 表明 综合 考虑 网 衣 和 和 鱼 类 的 作用 时 ,网 箱 下 游 流 速 衰 减 最 大 可 达 85% ,流速 约 在 6.5 倍 网 箱 直 
径 处 逐渐 恢复 ;与 仅 鱼 类 作用 和 仅 网 衣 作 用 相 比 ,网 箱 内 部 及 下 游 流速 衰减 程度 明显 增强 ,说 明 单 
二 < 独 只 考虑 鱼 或 网 部 不 能 准确 地 描述 鱼 类 养殖 网 箱 周 围 的 流 场 特性 , 须 将 二 者 结合 考虑 。 
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Characteristics of flow field around fish cages 
based on the porous media model 
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“Abstract :In order to explore the influence of fish cages on the flow field in the aquaculture area, based on 
the Reynolds time-average equation and Realizable k-e turbulence model, coupled with the porous media 
model ,a three-dimensional hydrodynamic model was established to study the flow field around fish cages. 
Three working conditions were carried out, namely, only the net effect, only the fish effect, and the com- 
bined effect of the net and the fish. The results show that when considering the effects of the net and the 
fish ,the attenuation of the flow velocity inside and downstream of the cage is significantly enhanced ,and 
the downstream flow velocity attenuation of the cage can reach a maximum of 85% ,and the starting posi- 
tion of the corresponding velocity recovery is about 6.5 times the diameter of the cage. Neither considering 
the fish or the net alone can accurately describe the characteristics of the flow field around fish cages,and 
the two factors must be considered together. The results of numerical simulation can provide generalized 
parameters for studies on the overall impact of aquaculture on the bay dynamics and further studies on the 
impact of aquaculture activities on water quality and environment. 
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近年 来 ,海水 网 箱 养殖 创造 了 巨大 的 经 济 效益 ; 
但 大 量 的 网 箱 聚 集 在 近 岸 海域 和 半 封 闭 海湾 "会 对 
养殖 区 的 水 流 起 到 阻碍 作用 ,降低 养殖 区 内 水 流 流 
速 ” ,从 而 引起 养殖 区 内 部 与 外 部 的 水 体 交 换 变 差 ， 
致使 养殖 区 内 水 环境 恶化 ,病害 频繁 发 生 , 养 殖 鱼 类 
大 量 死亡 ” 。 因 而 研究 养殖 网 箱 的 减 流 效应 有 利于 
网 箱 养殖 的 可 持续 发 展 。 

网 衣 系 统 是 网 箱 最 为 重要 的 结构 之 一 ,属于 柔 
性 小 尺度 \ 多 孔 结构 , 且 网 箱 中 的 养殖 鱼 类 数量 多 且 
具有 自主 游 动能 力 ,同时 受 养殖 活动 以 及 环境 的 影 
响 。 这 使 得 研究 鱼 类 养殖 网 箱 周围 流 场 的 影响 十 分 
困难 , 且 大 多 集中 在 现场 监测 和 物理 模型 试验 上 。 
王 辉 荣 " 通 过 对 湄 洲 湾 , 罗 源 湾 和 三 沙 湾 网 箱 养殖 
区 的 实测 数据 的 分 析 , 得 出 了 养殖 区 内 海水 基本 呈 
停 潇 状态 的 结论 。Johansson 等 中 在 挪威 的 4 个 大 
西 弹 刍 鱼 网 箱 养殖 区 进行 了 实地 监测 ,实测 数据 表 
明 网 箱 内 部 流速 降低 了 33% ~ 64% 。Gansel 等 
通过 物理 模型 试验 和 现场 试验 研究 了 空 网 箱 和 有 鱼 
网 町内 外 流 场 特性 ,结果 表明 相 比 于 空 网 箱 , 有 鱼网 
箱 会 在 网 箱 下 游 造 成 更 大 的 流速 衰减 。Plew 等 ” 
硬 帘 了 大 西洋 钙 鱼 对 圆 形 水 槽 中 水 流 的 影响 ,结果 
恬静 高 养殖 密度 下 水 槽 平均 流速 降低 了 57% 。 毕 
春 古 等 ”通过 水 槽 试验 研究 鲤鱼 幼 鱼 对 网 箱 周围 水 
动 小 特性 的 影响 ,结果 表明 对 于 两 个 串联 网 箱 ,来 流 
速度 越 大 , 鱼 群 对 网 箱 流速 的 影响 越 大 。 数 值 模拟 
方面 ,Tang 等 '” 利用 计算 流体 动力 学 模型 探究 了 养 
殖 忽 类 不 同 群 体 结构 游泳 速度 、 分 布 和 放养 密度 对 
网 箱 周 围 流动 动力 学 的 影响 。 由 于 网 箱 的 配 重 系统 
和 人 鲁 定 系统 起 到 固定 网 箱 的 作用 ,在 纯 水 流 作用 时 
可 忽略 网 衣 变 形 影 响 ,将 网 衣 看 作 一 个 整体 ,利用 多 
孔 介质 模型 进行 数值 模拟 。 目 前 已 有 较 多 的 研究 证 
明 多 孔 介质 模型 可 以 精确 地 模拟 单 片 网 衣 、 空 网 箱 
周围 流 场 特 性 '"” ,采用 多 孔 介质 模型 可 以 在 保证 
网 箱 尾 流 流 场 计 算 可 靠 性 的 基础 上 提高 计算 效率 。 

由 于 网 箱 中 的 养殖 鱼 数量 庞大 ,上 且 具 有 游泳 能 
力 和 独特 的 生活 习性 ,考虑 养殖 鱼 的 实时 分 布 特征 
及 游泳 状态 ,用 数值 模拟 的 方法 研究 养殖 鱼 对 网 箱 
周围 的 流 场 实时 影响 显然 是 不 现实 的 ,而 一 个 海水 
养殖 场 大 规模 的 养殖 行为 ,养殖 鱼 对 水 体 的 阻力 显 
然 不 能 忽略 。 从 统计 学 意义 上 ,可 以 认为 养殖 鱼 在 
特定 的 时 间 和 空间 是 均匀 分 布 的 。Gansel 等 ”通过 
实地 监测 发 现 网 箱 中 养殖 鱼 的 垂 向 分 布 有 一 定 的 规 
律 , 因 此 本 文 考虑 养殖 鱼 在 网 箱 中 的 垂 向 分 布 ,将 养 
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殖 鱼 覆盖 区 域 视 为 多 孔 介 质 区 域 ,利用 多 孔 介 质 模 
型 进行 模拟 计算 。 

本 文 基于 Ansys Fluent 平 台 , 视 网 箱 网 衣 以 及 鱼 
类 网 箱 内 部 为 多 孔 介 质 ,建立 了 人 研究 鱼 类 养殖 网 箱 
周围 流 场 的 三 维 水 动力 学 模型 。 在 验证 多 孔 介 质 模 
型 适用 性 的 基础 上 ,分 别 探究 了 网 箱 网 衣 、 鱼 类 以 及 
二 者 共同 作用 三 种 情况 下 网 箱 周围 的 流速 分 布 规 
律 。 单 体 鱼 类 养殖 网 箱 周 围 流 场 的 数值 模拟 可 为 网 
箱 养殖 海域 水 动力 模型 的 建立 提供 参数 支持 ,有 利 
于 进一步 评估 养殖 容量 以 及 优化 养殖 布局 。 


1 数学 模型 


1.1 控制 方程 
本 文 基 于 雷诺 时 均 方程 和 Realizable k-e 模型 ， 
耦合 多 孔 介 质 模型 对 网 箱 周围 流 场 进 行 求解 ,基本 

方程 为 
Wy + pu) = 


ot Ox, 9 a 
o(pu; o( puu,; 9 Ou, 
tp ) 9a(p 2 u pu )+s 
ot Ox; OX; Oxi\ Ax, , 
(2) 


式 中 :p 为 流体 密度 ;t 为 时 间 ;%; 为 空间 坐标 ;i 、w 
为 流体 速度 分 量 ;p 为 压强 ;ui \ 必 为 速度 分 量 的 脉动 
值 ;i,j =1,2,3 表示 坐标 分 量 ; 为 动力 黏度 ;$; 为 动 
量 方程 的 源 项 。 多 孔 介质 外 部 S; =0 ,多孔 介质 内 部 
动量 方程 的 源 项 为 

S = -~ (Dip + Cyz0 luly) (3) 


Dp, 0 0 Cc, 0 0 
Dj=|0 D 0lc=|o Cc, 0 
0 0 D 0 0 6， 


式 中 : 忆 , 为 法 向 黏 性 阻力 系数 ; 忆 为 切 向 黏 性 阻力 系 
数 ;C, 为 法 向 惯性 阻力 系数 ; C, 为 切 向 惯性 阻力 
系数 。 

由 于 网 衣 和 鱼 的 阻 水 作用 ,网 箱 前 后 存在 较 大 
的 压强 梯度 ,网 箱 附 近 的 水 流 运 动 变化 会 较为 激烈 ， 
Realizable fk-e 清流 模型 可 以 较 好 地 处 理 流速 变化 较 
大 的 流动 , 且 具 有 较 好 的 通用 性 "1; 因此 本 文选 用 
Realizable k-e 潮流 模型 进行 计算 。 
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1.2 多 孔 介 质 阻力 系数 计算 


1.2.1 和 鱼 类 阻力 系数 计算 

由 于 网 箱 中 的 鱼 数量 多 且 分 布 受 行为 习性 及 养 
殖 活动 影响 ,数值 模拟 困难 大 ;因此 将 网 箱 中 的 鱼 看 
作 一 个 整体 ,以 网 箱 中 的 体积 占 网 箱 相应 空间 体积 
的 比值 确定 多 孔 介质 模型 的 孔隙 率 , 将 网 箱 内 部 视 
为 多 孔 介质 区 域 ,实现 养殖 鱼 对 网 箱 周围 流 场 影响 ， 
如 图 1 所 示 。 


到 1 鱼 类 网 箱 简化 为 多 孔 介 质 模型 示意 图 

Fig.1 Sketch of porous media model 
秆 鱼 视 为 简单 均匀 的 多 孔 介 质 ,其 阻力 系数 可 
gun 公式 5 得 到 , 即 


2 
D, =D, 7 ea 
p & 
3.5(1—-e) (4) 
。 :a 
C Oe e3 


式 昌 : 忆 为 多 孔 介质 的 粒 径 ,由 养殖 鱼 的 体积 决定 ; 
eS -,e 为 孔隙 率 ,中 为 鱼 的 体积 分 数 ,由 网 箱 中 
的 全 的 体积 占 网 箱 相应 空间 体积 的 比值 确定 。 

生 2 ”网 衣 阻 力 系数 计算 

Patursson 等 "利用 8 片 网 衣 对 全 尺寸 单 体 网 
箱 赔 围 流 场 进行 数值 模拟 ,网 箱 下 游 流速 模拟 结果 
与 粟 测 数据 吻合 较 好 。 为 更 好 地 还 原 圆柱 形 网 箱 形 
状 特点 以 及 考虑 计算 效率 ,在 不 考虑 网 衣 的 变形 的 
情况 下 ,本文 将 网 箱 视 为 由 沿 圆周 均匀 分 布 的 12 片 
网 衣 组 成 ,通过 对 不 同 冲 角 下 各 个 平面 网 衣 模 型 多 
孔 介质 系数 的 设置 ,实现 对 单 体 网 箱 的 数值 模拟 。 
网 箱 简化 模型 如 图 2 所 示 。 


图 2 网 箱 简化 为 多 孔 介 质 模型 示意 图 
Fig.2 Sketch of porous media model 
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Patursson 等 '" 采用 最 小 二 乘法 ,根据 模型 试验 
结果 中 的 流速 与 阻力 的 关系 及 经 验 公 式 ,给 出 了 网 
衣 与 水 流 方向 垂直 时 ,多 孔 介质 系数 的 关系 式 为 

D,=1.11 x10’S’, +2.26x104S 
D,=3.97 x10°S2, +1.33 x10°5,, 
C, =70.257, +13.25,, 
C,=31.9S?, +6.935,, 
式 中 :5 为 网 衣 密 实 度 ;5,, =2d/a,d 为 网 线 直 径 ,a 
为 目 脚 长 度 。 

如 图 3 所 示 , 网 衣 与 水 流 成 一 定 冲 角 a。 此 时 ， 
可 将 多 孔 介 质 系数 按照 如 下 公式 进行 变换 "| 。 
_D,+D, D,-D 


(5) 


/ 1 
D’, 0 cos(2a’') 

C, + C C, C, , 
人 十 cos(2a’) 


(6) 


D,-D 
De 加 “sin(2a") 


C 一 C 
公用 “sin(2a') 


式 中 :ao =90° -a;D',、C',、D',、C' 为 变换 后 的 多 和 孔 
介质 系数 。 


上 


x 
pa < | 肢 长 度 a 
网 让 


图 3 网 衣 冲 角 示 意图 


Fig.3 Sketch of netting attack angle 


1.3 边界 条 件 及 数值 计算 方法 


数值 模拟 中 ,入 口 采用 速度 进口 边界 条 件 ,出 口 
采用 自由 出 流 边界 条 件 , 数 值 水槽 侧 壁 采用 对 称 边 
界 条 件 ,底面 采用 无 请 移 固 壁 边界 条 件 , 顶 面 采用 剪 
切 应 力 为 0 的 国 壁 边界 条 件 , 多 孔 介质 区 域 边界 条 
件 设置 为 多 孔 阶 跃 边界 条 件 。 采 用 有 限 体积 法 离散 
控制 方程 ,离散 格式 采用 二 阶 迎风 格式 ,压力 -速度 
耦合 选择 SIMPLE 算法 。 


1.4 模型 验证 


利用 Patursson “的 单 片 网 衣 的 试验 结果 对 网 
衣 模 型 进行 验证 ,如 图 4 所 示 。 试 验 水 覃 长 37 m , 宽 
3.66 m, 深 2.44 m; 试 验 网 衣 为 1 m x1 m 的 无 结 节 


尼龙 网 ,网 线 直径 为 2.8 mm, 目 脚 长 度 为 0.029 m。 
网 衣 由 钢 架 固定 在 水 覃 中央 ,位 于 水 面 以 下 0.73 m 
位 置 处 ,以 减少 网 衣 前 后 压力 差 产 生 的 表面 波 影 响 ; 
钢 架 可 绕 铬 垂 轴 旋 转 , 使 网 衣 与 水 流 方向 形成 不 同 
冲 角 。 数 值 模拟 的 计算 域 与 试验 水 槽 一 致 ,按照 网 
衣 冲 角 a 为 30" .45" .60" 和 90?* 四 种 工 况 进行 模拟 ， 
将 网 衣 视 为 多 孔 介 质 区 域 。 文 献 [8-9 ] 研究 表明 多 
孔 介 质 厚度 和 分 别 取 0.005 m.0.01 m 和 0.05m 时 
均 不 会 影响 网 衣 下 游 的 流速 衰减 , 故 本 文 多 孔 介 质 
区 域 厚度 取 0.05 m 进行 数值 模拟 。 经 计算 可 得 网 
衣 密 实 度 5,, 为 0.198。x、y 和 z 轴 分 别 代 表 流 向 、 
横向 和 垂 向 ,和信 口 为 均匀 来 流 , 速 度 为 U0, =0.5 m/s。 
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CS Fig.4 Sketch of model and measuring points 


己 提 取 不 同 冲 角 时 网 衣 中 心 后 2.5 m 处 的 流速 
U 避 计算 网 衣 的 阻力 系数 C, 和 升力 系数 C, ,模拟 
结果 与 试验 结果 对 比如 图 5 所 示 。 可 以 看 出 ,模拟 
结果 与 试验 结果 总 体 上 吻合 较 好 。 这 说 明 多 孔 介质 
模型 可 以 用 于 分 析 网 衣 以 及 网 箱 周围 的 流 场 特 性 。 
网 箱 中 的 鱼 .河道 中 的 挺 水 植被 等 都 会 对 水 流 
起 到 阻碍 作用 。 当 不 考虑 个 体 之 间 的 相互 作用 ,将 
其 看 成 一 个 整体 ,着 重 关注 尾 流 的 水 动力 特征 时 ,可 
将 网 箱 中 的 鱼 ,植被 等 整体 视 为 多 孔 介质 区 域 ,将 流 
动 区 域 中 这 些 固体 的 作用 近似 看 作 附加 在 流体 上 的 
分 布 阻力 进行 模拟 21 。 目 前 已 有 较 多 研究 表明 多 
孔 介质 模型 可 以 精确 地 模拟 植被 下 游 的 尾 流 特 
征 572] 。 本 文 利用 文献 [ 22 ] 的 挺 水 植被 试验 结 
对 养殖 鱼 类 多 孔 介 质 模型 进行 验证 。 如 图 6 所 示 ， 
试验 水 模 长 13 m, 宽 1.2m, 深 0.133 m;40 棵 挺 水 植 
物 交错 排列 形成 一 个 直径 也 为 0.22 m 的 圆 形 阵列 ， 
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植被 带 中 心 位 于 水 槽 中 轴线 上 ,高 度 与 水 槽 相同 , 单 
棵 植株 直径 a 为 0. 006 m。 本 文 计算 域 与 试验 水 权 
一 致 ,植被 区 域 视 为 多 孔 介 质 区 域 。 计 算得 到 植物 
固体 体积 分 数 为 0.03 ,入 口 流 速 为 Uo =0. 098 m/s。 
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Cj(C) 


(a) 流速 UU。 
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(b) CG, 和 GC, 


图 5 网 衣 的 阻力 系数 和 升力 系数 


Fig.53 Drag coefficient and lift coefficient of netting 
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(b) 侧 视图 
图 6 平面 计算 域 及 模型 示意 图 

Fig.6 Sketch of plane computing domain and model 
提取 沿 流向 方向 的 流速 U(y =0,z/H =0.5) 和 
沿 横向 方向 的 流速 V(y = D/2,z/H =0.5) 的 变化 曲 
线 与 试验 结果 进行 对 比 ,结果 如 图 7 所 示 。 可 以 看 
出 ,数值 模拟 结果 与 试验 值 较 为 吻合 ;流速 0 均 方 根 
误差 为 1.37% ,流速 了 均 方 根 误 差 为 0.79% 。 可 认 
为 多 孔 介 质 模 型 能 够 较 好 地 描述 多 孔 介 质 下 游 的 流 

场 结构 和 尾 流 特征 。 
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号 根据 海南 某 典型 养殖 场 实际 调查 资料 ,直径 D 
海地 ,高 度 妃 为 gm 的 圆 形 网 箱 , 养 鱼 数量 约 10 
万 是,6 月 份 单条 鱼 的 抽样 体重 均值 为 0.125 kg。 经 
计生 得 整个 网 箱 中 鱼 的 体积 分 数 由 为 0.22% 。 使 
用 酌 网 衣 网 线 直径 为 2.8 mm, 目 脚 长 度 为 0.029 m， 
经 针 算 得 网 衣 的 密实 度 ,为 0.2。 数 值 水 模 长 
519=m, 宽 95 m, 深 16 m。 网 箱 位 于 水 槽 中 央 , 网 箱 
中 已 距 离 速 度 人 口 38 m。 网 格 类 型 采用 非 结 构 化 四 
面体 网 格 ,将 网 衣 多 孔 介 质 区 域 按照 最 小 网 格 尺寸 
划分 网 格 ,再 将 其 周围 的 网 格 以 向 内 1. 1 倍 和 向 外 
1.2 倍 的 增长 比率 逐渐 增 大 ,同时 设置 最 大 网 格 尺 
寸 作为 网 格 尺寸 的 阔 值 。 通 过 网 格 无 关 性 检验 后 ， 
最 终 选 用 最 小 网 格 尺 寸 为 0. 025 m 的 网 格 进行 计 
算 , 局 部 网 格 划分 如 图 8 所 示 。 图 中 :xy 和 z 轴 分 
别 代表 流向 、 横 向 和 垂 向 ;红色 虚线 圈 表 示 网 箱 位 
置 。 入 口 为 来 流速 度 U, =0.5 m/s。 


2.1 仅 鱼 类 作用 时 网 箱 周围 流 场 特 性 研究 


当 养 殖 鱼 分 布 在 整个 网 箱 内 部 时 ,网 箱 周围 流 
速 分 布 如 图 9 所 示 。 图 中 黑色 点 画 线 区 域 表示 网 箱 
位 置 。 由 图 9 中 (a) 和 (c) 可 以 看 出 ,由 于 养殖 鱼 对 
水 流 的 阻碍 作用 ,网 箱 前 有 小 范围 的 流速 衰减 区 域 ; 
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网 箱 后 流速 衰减 区 域 范围 较 大 , 且 流速 骤 降 区 域 宽度 
逐渐 变 窄 ; 网 箱 两 侧 及 底部 水 流 流速 相应 增 大 。 这 是 
由 于 网 箱 中 的 鱼 的 阻 水 作用 ,水 流 会 在 网 箱 附 近 发 生 
绕 流 ,产生 横向 流动 和 垂 向 流动 ,从 而 网 箱 两 侧 及 底 
部 水 流 的 流向 发 生 明 显 偏转 。 这 与 Gansel 等 观测 
的 试验 现象 一 致 。 由 图 9 中 (b) 和 (d) 可 以 看 出 ,网 
箱 内 部 水 流 会 发 生 横向 和 垂 向 流动 ,产生 横向 的 流速 
V 和 垂 向 的 流速 WW; 横 向 流速 V 在 两 侧 达到 最 大 值 
0. 1 m/s, 为 来 流速 度 的 20% ; 垂 向 流速 玉 在 网 箱 
下 侧 达 到 最 大 ,为 0.075 m/s, 为 来 流速 度 U6 的 15%。 
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图 8 网 箱 局 部 网 格 示 意图 
Fig.8 Sketch of part meshing of net cage 
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(d) 垂直 截面 yD=0， 垂 向 流速 多 
图 9 网 箱 周围 流速 分 布 


Fig.9 Flow velocity distribution around the cage 


流向 流速 0 的 变化 曲线 如 图 10 所 示 。 可 以 看 
出 : 沿 x 方 向 ,网 箱 下 游 流速 先 减 小 ,后 逐渐 恢复 至 
稳定 值 ;由 网 箱 中 心 (y/D =0,z/ 有 = -0.5) 向 上 侧 
(y/D=0.25,z/H= -0.5;y/D =0.5,z/H= -0.5) 
靠近 时 ,流速 恢复 的 起 始 位 置 越 来 越 接近 网 箱 , 流 束 
衰减 程度 也 越 来 越 小 ; 穿 过 网 箱 中 心 沿 * 方向 流向 
流速 衰减 最 大 ,流速 衰减 可 达 80% , 约 在 8 倍 网 箱 直 
径 处 流速 开始 恢复 。 沿 方向 ,流速 呈现 网 箱 中 间 
低 ,两 侧 高 的 趋势 ;在 流向 流速 恢复 之 前 (x/D =5,z/ 
H= -0.5) ,横向 方向 流速 衰减 区 域 宽度 明显 大 于 
网 箱 直 径 ;在 流向 流速 沿 * 方向 恢复 之 后 (x/D = 10， 
15 ,20) , 尾 流 区 横 截面 上 流向 流速 会 在 网 箱 边缘 处 
恢复 成 来 流速 度 必 并 向 西 侧 逐渐 增 大 至 1. 09 们 来 
流速 度 。 沿 z 方向 , 网 箱底 部 流速 约 为 (1. 02 ~ 
LUD) ww , 比 来 流速 度 大 2% ~11% 。 
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图 10 流向 流速 UV 沿 不 同方 的 向 变化 趋势 
Fig. 10 Variation trend of flow velocity U 


along different directions 
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Gansel 等 '9 通 过 实地 监测 发 现 网 箱 中 养殖 鱼 类 
在 垂 向 分 布 上 有 明显 的 日 变化 现象 , 即 白天 鱼 群 多 
集中 在 网 箱 中 下 部 ,晚上 多 集中 在 网 箱 中 上 部 。 当 
鱼 群 集中 分 布 在 网 箱 中 上 部 或 者 中 下 部 时 ,中 = 
0. 33% 。 考 虑 养殖 鱼 类 在 垂 向 分 布 特点 对 流 场 的 影 
响 ,对 养殖 鱼 分 别 集中 分 布 在 网 箱 中 上 部 分 布 和 在 
网 箱 中 下 部 分 布 两 种 工 况 进行 数值 模拟 ,结果 如 图 
11 所 示 。 图 中 黑色 点 画 线 区 域 表示 网 箱 位 置 。 穿 
过 网 箱 中 心 提 取 沿 z 方向 流 癌 流速 U0 的 变化 曲线 如 
图 12 所 示 。 对 比 图 9(c) 和 图 11 可 以 看 出 , 当 养 殖 
鱼 类 集中 分 布 在 中 上 部 或 者 中 下 部 时 , 由 于 鱼 群 的 
体积 不 变 ,覆盖 区 域 变 小 ,因而 鱼 的 体积 分 数 会 相对 
增 大 ,给 水 流 造 成 更 大 的 阻力 ,水 流 的 垂 向 运动 加 
剧 ,在 网 箱 下 游 产 生 小 尺度 涡 旋 。 
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0 (a) 网 箱 中 上 部 
并 
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11 垂直 截面 网 箱 周围 流速 分 布 
Fig.11 Flow velocity distribution around the cage 
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图 12 流向 流速 U0 沿 z 方 向 变化 趋势 
Fig.12 Variation trend of flow velocity U along z direction 

由 图 12 可 以 看 出 : 鱼 群 集中 分 布 在 中 上 部 时 ， 

流速 U 在 网 箱 下 部 时 已 开始 恢复 ,并 在 网 箱底 部 时 

恢复 成 稳定 流速 1. 07 U6; 当 鱼 群 分 布 在 中 下 部 时 ， 

由 于 鱼 群 的 阻 水 作用 ,网 箱 上 部 水 流 流速 明显 增 大 ， 

可 达 1.08 U ,在 网 箱 下 游 流速 U 恢复 位 置 也 相对 于 

分 布 在 整个 网 箱 时 远离 网 箱 。 
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2.2 仅 网 衣 作 用 


图 13 为 仅 网 衣 作 用 时 网 箱 周围 的 流速 分 布 。 
图 中 黑色 点 画 线 区 域 表示 网 箱 位 置 。 由 图 13 中 (a) 
和 (b) 可 以 看 出 ,由 于 网 衣 的 阻 水 作用 ,同样 在 网 箱 
前 有 小 范围 的 流速 衰减 区 域 ;网 箱 后 流速 衰减 区 域 
范围 较 大 , 且 流 速 衰减 区 域 宽 度 大 致 与 网 箱 直径 相 
同 ;网 箱 两 侧 及 底部 水 流 流速 相应 增 大 ;由 流 线 可 以 
看 出 水 流 方向 始终 保持 与 人 口 流向 相同 ,说 明 网 衣 
的 阻 水 效应 并 不 会 引起 网 箱 下 游 及 底部 水 流 流 向 发 
生变 化 。 
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(c) 水 平 截面 z/H = -0.5 ,横向 流速 了 


(d) 垂 直 截 面 Y/D =0 , 垂 向 流速 丙 
图 13 网 箱 周围 流速 分 布 
Fig. 13 Flow velocity distribution around the cage 
由 图 13 中 (c) 和 (d) 可 以 看 出 ,网 衣 内 部 水 流 
会 发 生 横向 和 垂 向 流动 ,产生 横向 的 流速 了 和 垂 向 
的 流速 丈 ;横向 流速 了 在 两 侧 达到 最 大 值 0. 05 m/s， 
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为 来 流速 度 mm 的 10% ; 垂 向 流速 存在 网 箱 下 侧 达 
到 最 大 ,为 0.03 m/s, 为 来 流速 度 mm 的 6% 。 

流向 流速 UV 变化 曲线 如 图 14 所 示 。 可 以 看 出 : 
沿 x 方 向 ,网 箱 边 缘 处 (yAD =0.5,z/ 有 H= -0.5) 流 
速 0 变化 剧烈 , 尾 流 流速 减 小 后 迅速 开始 恢复 ; 穿 过 
网 箱 中 心 ,流速 先 减 小 至 0.78 倍 来 流速 度 , 并 于 2.5 
倍 网 箱 直径 处 流速 开始 恢复 。 沿 y 方向 ,网 箱 中 心 
流速 衰减 可 达 10% ;网 箱 尾 流 区 处 于 流速 衰减 区 域 
时 (x/D=1,z/ 有 = -0.5); 网 箱 内 部 尾 流 区 流速 先 
逐渐 恢复 ,后 趋 于 稳定 值 0.8U, , 而 网 箱 尾 流 区 处 于 
流速 恢复 区 域 时 (x/D = 10 ,20 ) ; 横 截 面 上 流速 呈现 
逐渐 减 小 后 趋 于 稳定 值 。 沿 z 方向 ,网 箱底 部 流速 
约 为 (1.01 ~1.03) Wo , 比 来 流速 度 大 1% ~3% 。 
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图 14 流向 流速 U 沿 不 同方 向 变化 趋势 
Fig. 14 Variation trend of flow velocity U 


along different directions 


2.3 和 鱼 类 和 网 衣 共 同 作用 


鱼 类 和 网 箱 共同 作用 时 ,养殖 网 箱 周围 流速 等 
线 图 如 图 15 所 示 。 图 中 黑色 点 画 线 区 域 表示 网 箱 
位 置 。 可 以 看 出 :其 流 场 特征 与 只 有 鱼 类 作用 时 相 
似 , 但 网 箱 内 部 和 下 游 流速 衰减 程度 明显 增强 ;网 箱 
内 部 水 流 的 横向 和 垂 向 流动 更 为 明显 ,横向 流速 V 
最 大 值 为 0. 14 m/s, 为 来 流速 度 0 的 28% , 垂 向 流 
速 WW 最 大 值 为 0.08 m/s, 为 来 流速 度 [的 16% 。 


wD 
| 
De 
© 
内 
© 
[> 
S 
ID 
Dn 


xID 


0 0.10 0.20 0.30 0.40 0.50 0.57m/s 


(a) 水 平 截面 zz 有 = -0.5, 流 向 流速 UV 


【=—) 


| 
[> 一 


So 
t 


So 
内 
So 
Ld 
Am 


xXxID 


0 0.10 0.20 0.30 0.40 0.50 0.57m/s 


De 
> 
LO 
守 
< 
© 
OO, 
CD 
© 
CD 
a 
© 
AN 
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(d) 垂 直 截 面 y/D =0, 垂 向 流速 多 
图 15 网 箱 周围 流速 分 布 
Fig.15 Flow velocity distribution around the cage 
取 穿 过 网 箱 中 心 流 向 流速 尽 变 化 曲线 与 相应 的 
只 有 鱼 作 用 和 只 有 网 衣 作 用 进行 对 比 , 结 果 如 图 16 
所 示 。 可 以 看 出 :考虑 鱼 和 网 共同 作用 时 流速 分 量 
U 变化 趋势 与 只 考虑 鱼 时 类 似 ,但 网 箱 内 部 和 下 游 
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流速 明显 比 仅 考虑 鱼 和 仅 考虑 网 时 低 , 网 箱 下 游 流 
速 衰减 可 达 85% ;相应 流速 恢复 的 起 始 位 置 会 靠近 
网 箱 , 约 为 6.5 倍 网 箱 直径 处 ; 相 比 于 只 有 鱼 类 作用 
时 ,网 箱 下 游 流速 增长 并 不 明显 ,但 流速 恢复 的 起 始 
位 置 略 远离 网 箱底 部 。 这 说 明 单独 只 考虑 域 或 网 都 
不 能 准确 的 描述 鱼 类 网 箱 周围 流 场 , 需 要 综合 考虑 
鱼 和 网 衣 的 作用 , 且 在 本 文 研究 的 密实 度 和 鱼 体积 
分 数 情况 下 , 鱼 对 水 的 减 流 作 用 明显 大 于 网 衣 。 


UU, 


xD 


UU, 


UU, 


ziH 
(c) :方向 


穿 过 网 箱 中 心 流 向 流速 UV 变化 趋势 
Variation trend of flow velocity changes 


across the center of the net cage 


本 文 基于 多 孔 介质 模型 ,建立 了 研究 鱼 类 养殖 
网 箱 周围 流 场 的 三 维 水 动力 模型 ,分 别 对 仪 网 衣 作 
用 、 仅 鱼 类 作用 和 二 者 共同 作用 3 种 情况 对 鱼 类 网 
箱 周围 水 动力 场 进行 了 模拟 ,分 析 了 网 箱 周 围 及 尾 
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流 特性 ,得 到 结论 如 下 : 

1) 多 孔 介质 模型 能 够 准确 地 模拟 大 尺度 鱼 类 养 
殖 网 箱 周围 的 流 场 特性 。 

2) 在 只 有 网 衣 作用 、 网 衣 密 实 度 5,。 =0.2 情况 
下 ,网 箱 周围 流 场 水 流 流向 不 发 生 明显 改变 ,横向 流 
速 了 和 垂 向 流速 WW 在 网 箱 周围 达到 最 大 ,最 大 值 分 
别 为 来 流速 度 的 10% 和 6% ; 网 箱 中 心 流速 衰减 可 
达 10% ,网 箱 下 游 流速 先 减 小 后 逐渐 增 大 ,最 大 流 
速 衰 减 约 22% ,并 于 2.5 倍 网 箱 直径 处 流速 开始 恢 
复 ; 网 箱底 部 水 流 流速 相应 增 大 , 约 比 来 流速 度 增 加 
1% ~3% 。 

3) 在 只 有 鱼 类 作用 、 鱼 的 体积 分 数 中 =0. 22% 

情况 下 ,由 于 网 箱 中 的 鱼 的 阻 水 效应 ,网 箱 周围 及 下 
游 水 流 流向 会 发 生 明显 偏转 , 且 网 箱 内 部 水 流 会 发 
生 弄 向 和 垂 向 运动 ,在 网 箱 周围 横向 流速 Y 和 垂 向 
流 妹 WW 会 达到 最 大 值 ,最 大 值 分 别 为 来 流速 度 的 
和 15% ;网 箱 下 游 流 向 流速 先 减 小 后 逐渐 恢复 
定 值 ,流速 衰减 最 大 可 达 80% ,并 约 在 8 倍 网 箱 
处 流速 开始 恢复 ;网 箱底 部 水 流 流速 相应 增 大 ， 
来 流速 度 增加 2% ~ 11% 。 考 虑 养殖 鱼 类 在 重 
厄 绍 布 特点 时 发 现 , 当 鱼 群集 中 分 布 在 网 箱 中 上 部 
或 部 中 下 部 时 ,网 箱 下 游 会 有 小 尺度 涡 环 产生 。 由 
而 群 的 体积 不 变 , 履 盖 区 域 变 小 , 因而 鱼 的 体积 分 
数 会 相对 增 大 ,给 水 流 造成 更 大 的 阻力 ,网 箱 下 游 流 
速 豆 减 更 剧烈。 
4 ) 在 鱼 和 网 衣 共 同 作用 时 , 流 场 分 布 特征 与 只 
盈 类 作用 时 类 似 。 但 经 过 对 比 可 以 看 出 ,网 箱 内 
下 游 流速 衰减 程度 明显 增强 ,网 箱 下 游 流速 衰 
减 最 大 可 达 85% ,相应 流速 恢复 的 起 始 位 置 会 靠近 
网 箱 , 约 为 6.5 倍 网 箱 直径 处 。 通 过 对 比 可 以 说 明 
单独 考虑 鱼 或 者 网 衣 作 用 均 不 能 准确 描述 鱼 类 网 箱 
周围 流 场 特性 ,需要 将 二 者 结合 考虑 。 
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